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Abstract 

We define and discuss a set of (4N - 4) parameters tliat can be used to analyse events 
in wliicli N jets liave been produced in fiigfi energy liadron-liadron collisions. Tliese 
multijet variables are the multijet mass and (4N - 5) independent dimensionless param- 
eters. To illustrate the use of the variables QCD predictions are presented for events 
with up to five jets produced at the Fermilab Tevatron Proton-Antiproton Collider. 
These QCD predictions are compared with the predictions of a model in which multi- 
jet events uniformly populate the N-body phase-space. 

Submitted to Physical Review D. 



1 



1 Introduction 



Large samples of events containing two or more jets have recently been recorded at the 
Fermilab Tevatron Proton- Antiproton Collider. Many of the observed events contain 
three-, four-, or even five-or-more jets |l|]. A comprehensive analysis of these multijet 
events would provide an interesting test of perturbative Quantum Chromo dynamics 
(QCD). In addition, a detailed understanding of the properties of multijet events pro- 
duced in high energy hadron-hadron collisions is important, firstly because multijet 
production is expected to be prolific in future high luminosity running at the Fermilab 
Proton-Antiproton Collider and at the Large Hadron Collider at CERN, and secondly 
because an understanding of QCD multijet production is required to facilitate the 
search for more exotic processes producing multijet events. For example, a detailed 
understanding of the properties of six-jet events at the Fermilab collider is likely to 
be important in the near future for the study of ti production and decay in the all 
hadronic channel. 

In the past, elegant analyses of two-jet and three-jet production have been published 
by the UAl ||, § and UA2 [g § collaborations at the CERN SppS CoUider and by the 
CDF [p, and DO collaborations at the Fermilab Tevatron Collider. There have 
also been analyses of events with more than three jets H, ^, |T0|. However, the analyses 
of events with four or more jets have not used a simple set of independent variables 
that (i) span the multijet parameter space, (ii) make it simple to interpret the observed 
event distributions within the framework of perturbative QCD, and (iii) make it is easy 
to compare the characteristics of events having N jets with the characteristics of events 
having for example (N-l-1) jets. In this paper we discuss a set of multijet parameters 
that satisfy these criteria. 

In choosing a set of multijet variables that span the multijet parameter space it 
should be noted that we can completely define a system of N particles in the N-body 
rest-frame by specifying 4N independent parameters, for example the 4N components 
of four-momentum. The N-body system would then be overspecified since momentum 
conservation provides us with three constraints. Furthermore, we can rotate the N-body 
system about the incoming beam direction without loosing any interesting information. 
Therefore, to describe the system we need only specify (4N - 4) parameters. We will 
take these parameters to be the N-body mass and (4N - 5) additional variables. We 
therefore introduce and discuss a set (4N - 5) dimensionless variables which, with the 
addition of the multijet mass, span the multijet parameter space. Our (4N - 5) multijet 
variables will provide a simple framework within which the properties of multijet events 
can be compared with QCD predictions. To illustrate the use of these variables we 
will compare predictions for the population of events in the multijet parameter space 
obtained from QCD matrix element calculations with the corresponding predictions 
from QCD parton shower Monte Carlo calculations, and from a model in which the 
events are uniformly distributed over the available phase-space. The QCD and phase- 
space calculations are described in section 2. In section 3 the analysis of two-jet events 
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is briefly discussed. The standard three-jet variables are reviewed and extended in 
section 4. Four-jet and five-jet variables are introduced and discussed in sections 5 and 
6. In section 7 the generalization of the multijet parameters to describe topologies with 
more than five jets is discussed. Finally, a summary is given in section 8. 

2 QCD and Phase-Space Predictions 

To illustrate the use of our multijet variables we will present and discuss various pre- 
dictions for the distribution of multijet events in the multijet parameter space. In 
particular we will consider two-jet, three-jet, four-jet, and five-jet events produced 
at the Fermilab Proton-Antiproton Collider operating at a center of mass energy of 
1.8 TeV, and compare predictions obtained from: (a) the HERWIG |0 QCD parton 
shower Monte Carlo program, (b) the NJETS [[12[ Leading Order (LO) QCD 2 ^ 
matrix element Monte Carlo program, and (c) a model in which events are distributed 
uniformly over the available N-body phase-space. 

2.1 Jet definitions and selection criteria 

The QCD and phase-space model predictions depend upon the algorithm used to define 
jets and the selection criteria used to define the data sample. To illustrate the use of our 
multijet variables we will take as an example jet definitions and event selection criteria 
recently used by the CDF collaboration to define a multijet data sample recorded at 
the Fermilab Proton-Antiproton Collider [0. Our predictions will therefore be for an 
existing data sample. Following the CDF prescription, jets are defined such that they 
satisfy the following: 

(i) Jet transverse energy Et > 20 GeV, where Et = E sin 6, E is the jet energy, and 

6 is the angle between the jet and the beam direction in the laboratory frame, 

(ii) I ?7 |< 3, where the jet pseudorapidity rj = — logtan(6'/2), and 

(iii) jet-jet separation AR > 0.9, where AR = [Arf + A</)^)^/^, and Ar] and Acf) are 
the differences in pseudorapidity and azimuthal angle between the two jets. 

With these jet definitions, the multijet event sample is defined by selecting events 
that satisfy the following: 

(a) Total transverse energy Et > 420 GeV, where the sum is over all jets with Et > 

20 GeV, 

(b) Multijet mass m > rrimin, and 

(c) the cosine of the leading-jet scattering angle cos6' < {cos 9)max where the leading- 

jet is defined as the highest energy jet in the multijet rest-frame. 
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Note that for two-jet events the J2 Et requirement selects events with jet Et > 
210 GeV. At fixed two-jet mass this results in an effective maximum allowed value 
of cos^. The values of nimin and {cos 9)max are chosen to restrict the parameter space 
to the region in which the YI^Et requirement is efficient. 



2.2 The HERWIG parton shower Monte Carlo calculation 



HERWIG [|ri|] is a leading-order QCD parton shower Monte Carlo program that in- 
cludes both initial- and final-state gluon radiation. HERWIG predictions can be 
thought of as LO QCD 2 — 2 predictions with gluon radiation, color coherence, 
hadronization, and an underlying event. We have used version 5.6 of the HERWIG 
Monte Carlo program together with a simple detector simulation that modifies the jet 
energies with a Gaussian resolution function: 

aE = 0.lE. (1) 

This is similar to the jet energy resolution function reported by the CDF collaboration 
In our HERWIG calculations we have used the CTEQIM structure functions 



and the scale = stu/2(s^+u^-|-t^). HERWIG generates 2 — > 2 processes above a 
specified p^"^*^ where p^"^'^ is the pt of the outgoing partons from the hard scatter before 
any radiation has occurred. We have set the minimum p^"^'^ to 60 GeV/c. Finally, the 
HERWIG Monte Carlo distributions discussed in this paper are inclusive. Hence, for 
a given jet multiplicity N, the generated events contribute to the distributions if they 
have at least N jets that pass the jet requirements. If there are more than N jets in a 
generated event, the multijet system is defined using the N highest Et jets. 



2.3 The NJETS QCD matrix element calculation 

The NJETS Monte Carlo program provides parton-level predictions based on the 
LO QCD 2 ^ N matrix elements. We have used the KMRSD- structure function 
parameterization's |T^ with the renormalization scale chosen to be the average pt 



of the outgoing partons. NJETS does not use a parton fragmentation model. Jet 
definitions and selection cuts are therefore applied to the final state partons. To enable 
a direct comparison between NJETS and HERWIG predictions we have smeared the 
final state parton energies in our NJETS calculations with the jet energy resolution 
function described above. 

In the following we will find that the NJETS and HERWIG predictions are generally 
in good agreement with one another. This suggests that the QCD predictions for the 
distributions discussed in this paper are probably not sensitive to reasonable variations 
in the choice of structure functions, scale, jet fragmentation (with the exception of 
the single-jet mass distributions), or underlying event. 
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2.4 Phase-Space model 



We have generated samples of Monte Carlo events for which the multijet systems uni- 
formly populate the N-body phase-space. These phase-space Monte Carlo events were 
generated with single-jet masses distributed according to the single-jet mass distri- 
bution predicted by the HERWIG Monte Carlo program. In addition, the multijet 
mass distributions were generated according to the corresponding distributions ob- 
tained from the HERWIG Monte Carlo calculation. Comparisons between the resulting 
phase-space model distributions and the corresponding HERWIG and NJETS Monte 
Carlo distributions help us to understand which multijet parameters are most sensitive 
to the behaviour of QCD multijet matrix elements. 



3 Two-Jet Variables 

We begin by briefly reviewing the variables that are often used in two-jet analyses 
0, 0, |[. Consider a system of two massless jets. The massless jet approximation is 
appropriate because at high center-of-mass energies single-jet masses are much smaller 
than two-jet masses (m2j). To describe a system of two massless jets in the two-jet 
rest-frame we need only two variables. In previous two-jet analyses these variables have 
often been chosen to be m2j and cos 6*, where 6* is the scattering angle between the 
incoming beam particles and the outgoing jets in the two-jet rest-frame. In defining 
cos 6** it must be remembered that in practice a two-jet system will always be produced 
together with a spectator system, and the incoming beam particles will not be coUinear 
in the two-body rest-frame. Hence, following the convention of Collins and Soper |15 



6* is taken to be the angle between the outgoing jets and the average beam direction. 
Consider the process 1 + 2 — 3 + 4. The center-of-mass scattering angle is defined: 

n* — AV ■ ^3 fc.\ 

I -TyU/ II -TS I 

where 

'Pav = i^i-i^2, (3) 

and we define particle 1 as the incoming interacting parton with the highest energy in 
the laboratory frame. 

NJETS and HERWIG QCD Monte Carlo predictions for the maj and cosO* dis- 
tributions are shown in Figs. |l] and ^ respectively for two-jet events produced at the 
Fermilab Proton- Antiproton Collider satisfying the requirements m2j > 550 GeV/c^ 
and I cos 6'* | < 0.6. Note that in the HERWIG Monte Carlo calculation the jets 
acquire mass in the fragmentation process, whereas in the NJETS calculation jets are 
identified with massless partons. Hence the agreement between the HERWIG and 
NJETS predictions reflects the validity of the massless jet approximation. The pre- 
dicted cos^* distributions are similar to the angular distribution expected at LO for 
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qq qq scattering which is not very different from the well known Rutherford 
scattering form: 

'^'^ i (l-cosr)-2. (4) 



d cos 9* 



Hence, the cos 6'* variable has some nice features. Firstly, the LO QCD prediction 
for the cos 6'* distribution is well known and is similar, although not identical, to the 
Rutherford scattering distribution. Secondly, the phase-space density is independent of 
cos^*. Therefore the measured cos 6* distribution depends upon the underlying 2 —>■ 2 
matrix element in a very direct way. 

To prepare for the analysis of systems with many jets in the final state it is useful 
to extend the two-jet variables to describe two-jet systems with massive final state 
jets. To do this we must specify two additional parameters. Obvious choices are the 
final state single-jet masses and m^. We prefer to use dimensionless variables, and 
therefore define: 



and 



h - ^, (5) 

m2j 



fi = , (6) 

m2j 



where the jets are ordered in the two-body rest-frame such that > E4, and hence 
/s > /4. The HERWIG predictions for the /s and distributions are shown in Fig. |^. 
Note that /s and /4 are typically of order 0.05 to 0.1, and hence single-jet masses can 
be neglected for many purposes. 

We conclude by noting that we have defined four variables that specify a two-jet 
system in the two-body rest-frame: m2j, cos^*, /s, and f^. 



4 Three-Jet Variables 

In the standard three-jet analysis used by the UAl collaboration 0, and later by the 
CDF 1^ and DO collaborations, five variables are chosen that specify the system 
of 3 massless particles in the three-body rest-frame. The first of these variables is the 
three-jet mass {rrisj). The NJETS and HERWIG predictions for the rrisj distribution 
are shown in Fig. ^ to be in good agreement with each other. The predicted m^j 
distributions have also recently been shown to be in good agreement with the observed 
CDF msj distribution To complete the description of the three-jet system four 
additional dimensionless variables are defined that, together with 777.3 j, span the three- 
body parameter space. In defining the three-jet parameters it is traditional to label 
the outgoing jets 3, 4, and 5, and order the jets such that E3 > E4 > E^, where Ej 
is the energy of jet j in the three-jet rest-frame. The traditional three-jet variables 
employed are X3, X4, cos 6*3, and ips, which are defined: 
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(i) X3, the leading-jet energy fraction, normalized: 



X = ^ ^3 ^ 2_E3 

(ii) X4, the next-to-leading jet energy fraction, normalized: 

(iii) cos 6'3, defined in the three-jet rest-frame as the cosine of the leading-jet scattering 
angle (see Fig. |^) : 

a — AV ■ P3 ,f.s 

I -T Ay II -T 3 I 

(iv) tps, defined in the three-jet rest-frame as the angle between the three-jet plane 
and the plane containing jet 3 (the leading jet) and the average beam direction 
(see Fig. H) : 

_ (7^3 X Vav) ■ (P4 X 7^5) .... 

= 17^ x7^ 117^ x7^ I • ^ ^ 

I -r 3 X r'AV II -r 4 X -Ts | 

Predictions for the X3-, X4-, cos6'3-, and Vs^distributions are shown in Fig. ^ 
for three-jet events produced at the Fermilab Proton-Antiproton Collider that satisfy 
the requirements rrisj > 600 GeV/c^, | cos 6^3 |< 0.6, and X3 < 0.9. These selection 
criteria are used to restrict the parameter space to the region for which the Y^Et 
requirement is efficient and to ensure that the jets in the three-jet sample are well 
measured. The first and second three-jet parameters (X3 and X4) are Dalitz variables, 
normalized so that X3 + X4 -|- X5 = 2. Momentum conservation restricts the ranges of 
the Dalitz variables (for massless jets 2/3 < X3 < 1 and 1/2 < X4 < 1). The phase- 
space density is uniform over the kinematically allowed region of the (X3, X4)-plane, 
and hence the phase-space model predictions for the X3- and X4-distributions can be 
easily understood. Note that the QCD predictions for the X3- and X4-distributions 
are similar to those of the phase-space model. We might have expected the QCD 
calculations to predict an enhanced event rate as X3 — >■ 1 and the three-jet system 
therefore approaches a two-jet configuration. However, in practice the algorithm used 
to define jets and the experimental requirements used to select well measured three-jet 
events restrict the measured three-jet topologies to those that populate regions of the 
three-body phase-space where the matrix element varies only slowly over the (X3, X4)- 
plane. The third and fourth three-jet parameters (cos6'3 and ip^) are angular variables. 
The phase-space density is uniform in cos 6'3-space, ?/'3-space, and is also uniform in 
the (cos 6*3, '?/'3)-plane. Indeed, the phase-space model does predict a uniform cos ^3 
distribution. The phase-space model prediction for the distribution is not quite 



7 



uniform, there being a slight depletion of events as t/^s — or vr. This depletion is 
primarily a consequence of the minimum Et requirement used to define jets. We would 
expect the QCD predictions for the two angular distributions to be very different from 
the phase-space model predictions. In particular we might expect that the leading- 
jet angular distribution would be similar, although not identical, to the LO qq qq 
scattering form. Indeed, this is seen to be the case for both the NJETS and HERWIG 
QCD predictions (Fig. ^c). We might also expect the initial-state radiation pole in the 
QCD matrix element to result in an enhanced rate of three-jet events for topologies in 
which the angle between the beam direction and the three-jet plane is small. Hence, 
we would expect the ips distribution to be peaked towards and vr. This is also evident 
in the HERWIG and NJETS predictions. 

To prepare for the analysis of events with more than three jets we now wish to 
extend the three-jet variables to describe a system of three massive particles in the 
three-body rest-frame. To do this we must specify an additional three parameters, 
which we take to be the single-jet mass fractions f^, f^, and /s defined: 

(a) /s, the leading-jet mass divided by the three-jet mass: 

/, - ^, (11) 

(b) /4, the next-to-leading-jet mass divided by the three-jet mass: 

, _ m4 
h = ) 

(c) /s, the third-to-leading-jet mass divided by the three-jet mass: 



(12) 



/= - ^. (13) 

HERWIG predictions for /s, /4, and /s are shown in Fig. ^ Note that fj is typi- 
cally less than or of order 0.1, and hence single-jet masses can be neglected for many 
purposes. 

We conclude by noting that we have defined eight variables that specify a three-jet 
system in the three-body rest-frame: msj, X3, X4, cos 6*3, t/'s, /s, /4, and /s. 



5 Four-Jet Variables 

To completely describe a system of four jets in the four-body rest-frame we must specify 
twelve independent parameters. We will choose the four-jet mass (m4j) and eleven 
dimensionless variables that span the four-body parameter space. We have chosen a 
set of four-jet variables that, for four-jet configurations that approach a three-body 
topology, reduce to the three-jet variables discussed in the previous section. This will 
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make it possible to compare the characteristics of four-jet events with the corresponding 
characteristics of three-jet events. 

The four-jet variables are shown schematically in Fig. p. We begin by reducing 
the four-jet system to a three-body system by combining the two jets with the lowest 
two-jet mass. We will label the two jets we combine A and B with Ea > Eb, where 
Ea and Eb are the jet energies in the four-jet rest-frame. The resulting three-body 
system can be completely specified using our three-jet variables: Xy, X4/, cos 6^3/, ipy, 
/s') /4') and /s/. Note that we order the three bodies in the three-body rest-frame so 
that E3/ > E^i > E^', and use a nomenclature in which primed labels denote objects 
after two jets have been combined. Hence one of the three primed objects will be the 
two-jet system (AB). Explicitly, X3/, X4/, cos 6*3/, ip^i, fy, f^i, and f^' are defined: 

(i) X3/, the fraction of the three-body energy taken by the leading object, normalized: 

2 E'li 2 E-v 

X3/ = — — = , 14 

E-y + Ei, + E^. m4j ' ^ ^ 

(ii) X4/, the fraction of the three-body energy taken by the next-to-leading object, 
normalized: 

X4, = ^ ^ ^ ^ ,^15) 

Ey + E^i + E^i m^j 

(iii) cos 6*3/, the cosine of the leading-body scattering angle: 

cos ^3' = , (16) 

I -r Ay II "3' I 

(iv) ipy, the angle between the three-body plane and the plane containing object 3' 
(the leading body) and the average beam direction: 

_ (:^3' X 'Pay) ■ (Pa' X ?50 
cosips' = =^ ^ ^ — , (17) 

I P3' X Pay II P4' X P5. | 

(v) /s', the mass of the leading object divided by the four-jet mass: 

/. = ^. (18) 

(vi) /4/, the mass of the next-to-leading object divided by the four-jet mass: 

U ^ ^, (19) 

(vii) /s/, the mass of the third-to-leading object divided by the four-jet mass: 

/. - ^. (20) 



9 



The NJETS and HERWIG predictions for the m^j distribution are shown in Fig. ^ 
for four-jet events produced at the Fermilab Proton- Antiproton Colhder satisfying the 
requirements m4j > 650 GeV/c^, | cos 631 \< 0.8, and X31 < 0.9. The QCD predictions 
for the X3/-, X4/-, cos 6^3/-, and t/'s/ -distributions are compared with the phase-space 



model predictions in Fig. [1^. There is reasonable agreement between the HERWIG 
and NJETS predictions for all of these distributions. The QCD predictions for the 
Xy- and X4/-distributions are not very different from the predictions of the phase- 
space model. In contrast, the NJETS and HERWIG cos6'3/- and ?/'3/-distributions are 
very different from the more uniform phase-space model predictions. It is interesting 
to compare these distributions with the equivalent distributions for three-jet events 
(Fig. ^). The QCD and phase-space model predictions for the four-jet distributions 
are similar but not identical to the corresponding distributions for three-jet events. 
Note that (1) in comparing the phase-space model predictions for the X3- and X3/- 
distributions we see that the predicted Xy distribution is depleted at large Xy, and 
(2) in comparing the phase-space model predictions for the X4- and X4/-distributions 
we see that the predicted X^/ distribution is distorted at large X4/. These differences 
can be qualitatively understood by noting that if 4' or 5' is the (AB)-system and hence 
massive then Xy < 1 even if 4' and 5' are coUinear. It should also be noted that the 
phase-space model cos 6*3' distribution is slightly depleted at small | cos 6*3' | and the 

distribution is slightly depleted for values of ip^' close to and vr. These features 
are consequences of the minimum jet-jet separation requirement AR > 0.9, and the 
minimum jet transverse energy requirement Et > 20 GeV. 

The HERWIG predictions for the normalized single-jet masses /-,/ are shown in 
Fig. 0. They exhibit peaks close to /,/ = 0.05 which reflect the finite single-jet masses 
resulting from the HERWIG fragmentation model, and long tails at larger values of fj' 
which reflect the contributions from the combined (AB)-systems. Note that although 
single jets are massless in the NJETS calculation, the NJETS program does predict 
the contribution to the /^'-distributions from the combined (AB)-systems, and indeed 
the NJETS and HERWIG predictions are in good agreement at large fj'. 

To complete our description of the four-jet system we must now specify four addi- 
tional parameters that describe the two-jet (AB)-system. To describe the (AB)-system 
we choose: 

(a) Ja, the mass of jet A divided by the four-jet mass: 

h ^ ^, (21) 



(b) fB, the mass of jet B divided by the four-jet mass: 



fs ^ ^, (22) 
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(c) Xa, defined in the four-jet rest-frame as the fraction of the energy of the (AB)- 
system taken by the leading jet: 

(d) ipAB^ defined in the four-jet rest-frame as the angle between (i) the plane contain- 
ing the (y4i?)-system and the average beam direction, and (ii) the plane containing 
A and B (see Fig. |^). The prime reminds us that in order to define ip'^B we have 
combined two jets to obtain the (AB)-system. Note that: 

. X (|.,x ?,,,,) ^ 

\ Pax Pb II P ab x P AV I 

The predicted /a- and /s- distributions are shown in Figs. |T2] (a) and [T^ (b) respec- 
tively. The typical values of and predicted by the HERWIG fragmentation model 
are less than or of order 0.05. The predicted Xa distributions are shown in Fig. |l^ (c). 
The NJETS and HERWIG QCD calculations yield harder Xa distributions than the 
corresponding distribution predicted by the phase-space model. Presumably this re- 
flects the presence of the soft gluon radiation pole in the QCD matrix element. To gain 
some insight into the shape of the phase-space model prediction for the Xa distribution 
consider a system of four massless particles labelled randomly i, j, k, and 1. If we define 
Xi = Ei/{Ei + Ej), then the phase-space prediction for the distribution of events as a 
function of Xj is given by: 

dN_ ^ J ^25) 

dXi Xf Xf 

This function is already quite similar to the phase-space model prediction shown in 
Fig. |12| (c), which is obtained by requiring that the (AB)-system is the lowest mass pair, 
and taking account of finite single-jet masses and experimental selection requirements. 
Finally, the predicted iP'ab distributions are shown in Fig. |12| (d). The NJETS and 
HERWIG predictions for the iP'ab distribution are in agreement with one another. The 
slight decrease in the population of events predicted by the phase-space model as iP'ab 
approaches or vr is a consequence of the minimum jet Et requirement. 

We conclude by noting that we have defined twelve variables that specify a four-jet 
system in the four-body rest-frame: m4j, X3/, X4/, cos 6*3/, ^3/, /3/, fi>, /g/, /a, /b, Xa, 
and Vab- 

6 Five-Jet Variables 

To completely describe a system of five jets in the five-body rest-frame we must specify 
sixteen independent parameters. We will choose the five-jet mass (m^j) and fifteen 



11 



dimensionless variables that span the five-body parameter space. We have chosen a set 
of five-jet variables that, for five-body configurations that approach a four-body topol- 
ogy, reduce to the four-jet variables discussed in the previous section. Furthermore, for 
five-body configurations that approach a three-body topology, our five-jet parameters 
reduce to the three-jet variables discussed previously. Thus we will be able to compare 
the characteristics of five-jet events with the corresponding characteristics of three-jet 
and four- jet events. 



The five-jet variables are shown schematically in Fig. |I3[ We begin by reducing 
the five-jet system to a four-body system by combining the two jets with the lowest 
two-jet mass. We will label the two jets we combine C and D, with Ec > Eo, where 
Ec and Ej:, are the jet energies in the five-jet rest-frame. We can then further reduce 
the resulting four-body system to a three-body system by combining the two bodies 
with the lowest two-body mass. We will label the two objects we combine A' and B', 
with Ea' > Eb'- The resulting three-body system can be completely specified using 
our three-jet variables: X3", X4//, cos 9yi, ipyi, fyi, f^ir, and /s". Note that we order 
the three bodies such that Eyi > E^n > E^n, and use a nomenclature in which doubly 
primed labels denote objects after two operations in which the two bodies with the 
lowest two-body mass have been combined. One of the three doubly primed objects 
will be the (A'i?')-system. Explicitly, Xs», X4//, cos 6^3//, ips", fs", f^", and /s" are 
defined: 

(i) Xyt, the fraction of the three-body energy taken by the leading body, normalized: 

2 E'iii lEvi 
E^ii + E^ii + E^n m^j 

(ii) X4/', the fraction of the three-body energy taken by the next-to-leading body, 
normalized: 

V - 2 E^n 2-E4// 

= p , p , p - — — ' \^') 

(iii) cos6'3'/, the cosine of the leading-body scattering angle: 

cos ^3" = — (28) 

I ^ AV \ \ "3" I 

(iv) ipyi, the angle between the three-body plane and the plane containing object 3" 
(the leading body) and the average beam direction: 

= ^7? 117^ x7^ I ' ^ ^ 

I J 3" X r AV II -14" X J 5" I 

(v) fy, the normalized mass of the leading object: 



h" - ^, (30) 

m5j 
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(vi) f^n, the normalized mass of the next-to- leading object: 



- ^, (31) 

(vii) /s//, the normalized mass of the third-to-leading object: 

h" - (32) 

The NJETS and HERWIG predictions for the m^j distribution are shown in Fig. 
for five-jet events produced at the Fermilab Proton- Antiproton Collider and satisfying 
the requirement m^j > 750 GeV/c^. The QCD predictions for the X3//-, X4//-, cos 6^3//-, 



and ?/'3//-distributions are compared with the phase-space model predictions in Fig. [15. 
The predicted distributions are qualitatively similar to the equivalent four-jet distri- 
butions shown in Fig. |l^. Note that the QCD predictions for the cos distribution 
are remarkably similar to the simple LO qq — > qq angular distribution. 

The HERWIG predictions for the normalized single-jet masses fj" are shown in 
Fig. O. Once again, the HERWIG and NJETS distributions are in agreement at large 



mass fractions. 

We must now specify the intermediate four-body system. In analogy with the four- 
jet analysis we will do this by specifying four additional dimensionless variables that 
describe the (y4'i?')-system. We choose /a', fs', ^a', and V'a'B'' defined: 



(a) /a', the normalized mass of object A 



I. 



h. - ^, (33) 

(b) fs'i the normalized mass of object B': 

/b. - ^, (34) 

(c) Xa'i defined in the five-jet rest-frame as the fraction of the energy of the [A'B')- 
system taken by the leading body: 

Xa' = ^ , (35) 

Ea' + Eb' 

(d) i'A'B'y defined in the five-jet rest-frame as the angle between (i) the plane con- 
taining the (74'i?')-system and the average beam direction, and (ii) the plane 
containing A' and B' (see Fig. |13D. Note that: 

,// C^A' X P B') ■ [P A'B' X P Av) tOR\ 

cos^A,^, = ^ — ^ . (36) 

\ P A' ^ P B' W P A'B' ^ P AV \ 
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The predicted distributions of these variables are shown in Fig. |T7[ The HERWIG 
predictions for the /a'- and /^'-distributions peak at values of about 0.02 and have 
long tails associated with composite A' or B' systems. The tails are accounted for by 
the NJETS predictions. It is interesting to compare the Xa'- and ?/'^,^/-distributions 
with the corresponding four-jet distributions (Figs. ^ (c) and[l^ (d) respectively). The 
QCD and phase-space model predictions for the five-jet distributions are qualitatively 
similar to the corresponding four-jet distributions. Note that the HERWIG and NJETS 
predictions are in agreement with one another. 

Finally, to complete our specification of the five-jet system we must define a further 
four variables that describe the two-body (CD)-system. We choose fc, fn, Xc, and 
ipcjy, defined: 



(a) fc, the normalized mass of jet C: 



fc ^ ^. (37) 

(b) fu, the normalized mass of jet D: 

fo ^ ^, (38) 

(c) Xc, defined in the five-jet rest-frame as the fraction of the energy of the (CD)- 
system taken by the leading jet: 

Xc - TT^, (39) 

(d) ipcD^ defined in the five-jet rest-frame as the angle between (i) the plane contain- 
ing the (C-D)-system and the average beam direction and (ii) the plane containing 
C and D (see Fig. |1|). Note that: 



COS ^CD = -j^ ^ -■ (40) 

\ Pc'X Pd II PcD X P AV 



The predicted distributions of these variables are shown in Fig. |T8[ The HERWIG 
predictions for the /a'- and /^'-distributions peak at values less than 0.02. Note 
that the QCD predictions for the Xc distribution are harder than the corresponding 
phase-space model prediction, whilst the QCD predictions for the ^/'^^-distribution 
are similar to the corresponding phase-space model prediction. 

We conclude by noting that we have defined sixteen variables that specify a five-jet 
system in the five-body rest-frame: msj, X3//, X4//, cos 6*3", ■j/'s", /a", /4", /s", fA', fs', 
Xa', i^A'B'^ fc, Id, Xc, and 
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7 Generalization to Events with Six or More Jets 



A list of the multijet variables described in the preceeding sections is given in Table |l|. 
The extension of the variables to describe multijet systems with more than five jets is 
straight forward. As an example the variables required to describe a six-jet event are 
also listed in Table |I|. In general, to describe an event containing N jets we use the mass 
of the N-jet system plus (4N-5) dimensionless variables. To define the dimensionless 
variables we proceed by reducing the N-jet system to a three-body system. This is 
done in (N-3) steps. In each step the two bodies with the lowest two-body mass are 
combined by adding the two four- vectors. The resulting three-body system is described 
by specifying seven parameters, namely the normalized masses of the three bodies (e.g. 
/s, /4, and /s), the Dalitz variables for the two leading bodies (e.g. X3 and X4), the 
cosine of the leading-body scattering angle (e.g. cos ^3), and the angle between the 
three-body plane and the beam direction (e.g. ip^). To complete the description of the 
N-jet system we must then specify an additional four parameters for each step in which 
two bodies were combined. These parameters are the normalized masses of the two 
bodies (e.g. and /b), the fraction of the two-body energy taken by the leading body 
(e.g. Xa)i and the angle defined in the N-jet rest- frame between the plane containing 
the two-body system and the beam direction and the plane defined by the two bodies 
(e-g V'ab)- 

8 Summary 

We have defined a set of (4N - 4) parameters that can be used to analyse events in 
which N jets have been produced in high energy hadron-hadron collisions. These mul- 
tijet parameters (i) span the multijet parameter space, (ii) facilitate the interpretation 
of observed event distributions within the framework of perturbative QCD, and (iii) 
make it is possible to compare the characteristics of events having N jets with the 
characteristics of events having for example (N-l-1) jets. 

To illustrate the use of the multijet variables described in this paper we have dis- 
cussed QCD and phase-space model predictions for three-jet, four-jet, and five-jet 
events produced at the Fermilab Proton-Antiproton Collider. For this particular ex- 
ample we note that, apart from small effects that can be ascribed to the absence of a 
fragmentation model in the NJETS calculation, the complete LO QCD matrix element 
predictions for each of the multijet parameter distributions discussed in the preceed- 
ing sections are well described by the parton shower Monte Carlo calculation. Thus 
it appears that even when there are five hard partons in the final state a good ap- 
proximation to the LO QCD matrix element is given by 2 — >^ 2 scattering plus gluon 
radiation. This is of interest because the complete LO matrix element calculation is 
not at present available for topologies with more than five final state partons. Hence 
for the analysis of events with six or more jets we must rely on parton shower Monte 
Carlo calculations, or on other approximations to the QCD matrix element. 
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Finally, the multijet variables discussed in this paper have been selected to em- 
phasize simple to interpret quantities (masses, energy fractions, scattering angles, and 
planarity-type angles) . Wherever possible we have tried to select parameters for which 
the phase-space model distributions are simple to understand. Experimental require- 
ments used to select well measured N-jet events necessarily distort some of the pre- 
dicted distributions. However, for the example discussed in this paper, we note that 
for most parameters the experimental selection criteria result in modifications to the 
phase-space model predictions that are modest and are limited to small regions of the 
parameter space. The observed N-jet distributions should therefore directly reflect the 
dynamics of the underlying multijet matrix element. 

We are grateful to Walter Giele for many interesting discussions. This work was 
supported by the U.S. Department of Energy, and the Ministry of Science, Culture and 
Education of Japan. 
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Table 1: Summary of the (4N-4) multijet variables for N = 2, 3, 4, 5, and 6. 
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CM 




m,, (GeV/c') 

Figure 1: Predicted two-jet mass distributions for two-jet events produced at the Fermi- 
lab Proton-Antiproton Collider. HERWIG (points) compared with N JETS (histogram) 
after applying the requirements of m2j > 550 GeV/c^ and | cos 6** |< 0.6. 




Icosi^l 

Figure 2: Predicted | cos^* | distributions for two-jet events produced at the Fer- 
milab Proton-Antiproton Collider that satisfy the requirements m2j > 550 GeV/c^ 
and I cos^* |< 0.6. The HERWIG prediction (points) is compared with the NJETS 
prediction (histogram), and the LO QCD prediction for qq — > qq scattering (curve). 
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Figure 3: The HERWIG Monte Carlo predictions for the distributions of leading and 
next-to-leading single-jet-mass fractions for jets in two-jet events produced at the Fer- 
milab Proton- Antiproton Collider that satisfy the requirements m2j > 550 GeV/c^ 
and I cos^* |< 0.6. 
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Figure 4: Predicted three-jet mass distributions for events produced at the Fermilab 
Proton- Antiproton Colhder that satisfy the requirements rrisj > 600 GeV/c^, X3 < 
0.9, and | cos ^3 |< 0.6. HERWIG predictions (points) are compared with NJETS 
predictions (histogram) . 



THREE-BODY REST FRAME 




Figure 5: Schematic definition of angles used to describe the three-jet system in the 
three-jet rest-frame. 
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Figure 6: Predicted distributions of the three-jet variables defined in the text for 
three-jet events produced at the Fermilab Proton- Antiproton Colhder that satisfy the 
requirements rrisj > 600 GeV/c^, < 0.9, and | cos6'3 |< 0.6. HERWIG predictions 
(points) are compared with NJETS predictions (histograms) and the phase-space model 
predictions (solid curves) for (a) X3, (b) X4, (c) cos ^3, and (d) ■03- The broken curve 
in the cos 63 figure is the LO QCD prediction for qq qq scattering. 
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Figure 7: HERWIG Monte Carlo predictions for the singlc-jct mass-fraction distribu- 
tions for jets in three-jet events produced at the Fermilab Proton- Antiproton CoUider 
that satisfy the requirements rrisj > 600 GeV/c^, < 0.9, and | cos6'3 |< 0.6. 
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Figure 8: Schematic definition of angles used to describe the four-jet system in the 
four-jet rest-frame. 
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Figure 9: Predicted four-jet mass distributions for events produced at the Fermilab 
Proton- Antiproton Collider that satisfy the requirements > 650 GeV/c^, X3/ < 
0.9, and | cos ^3/ |< 0.8. HERWIG predictions (points) are compared with NJETS 
predictions (histogram) . 
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Figure 10: Predicted distributions of three-body variables described in the text for 
four-jet events produced at the Fermilab Proton- Antiproton CoUider that satisfy the 
requirements rriij > 650 GeV/c^, X3/ < 0.9, and | cos% |< 0.8. The HERWIG 
predictions (points) are compared with NJETS predictions (histograms), and with the 
phase-space model predictions (solid curves) for (a) X^/, (b) X4/, (c) cos 6*3/, and (d) 
■03/. The broken curve in the cos 6*3/ figure is the LO QCD prediction for qq qq 
scattering. 
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Figure 11: The predicted distributions of single-jet mass-fractions for jets in four-jet 
events produced at the Fermilab Proton- Antiproton Colhder that satisfy the require- 
ments m4j > 650 GeV/c^, X^, < 0.9, and | eos^g/ |< 0.8. HERWIG predictions 
(points) are compared with NJETS predictions (histograms). 
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Figure 12: The predicted distributions of the four-jet variables describing the (AB)- 
system for four-jet events produced at the Fermilab Proton- Antiproton CoUider that 
satisfy the requirements m4j > 650 GeV/c^, < 0.9, and | cos ^3/ |< 0.8. The 
HERWIG predictions (points) arc compared with NJETS predictions (histograms), 
and the phase-space predictions (curves) for (a) Ja, (b) Jb, (c) Xa, and (d) iP'ab- 
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Figure 13: Schematic definition of angles used to describe the five-jet system in the 
five-jet rest-frame. 
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Figure 14: Predicted five-jet mass distributions for five-jet events produced at tlie 
Fermilab Proton- Antiproton Collider. HERWIG predictions (points) compared with 
NJETS predictions (histogram). 
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Figure 15: Predicted distributions of three-body variables for five-jet events pro- 
duced at the Fermilab Proton- Antiproton Colhder that satisfy the requirement m^j > 
750 GeV/c^. HERWIG predictions (points) are compared with NJETS predictions 
(histograms) and the phase-space predictions (sohd curves) for (a) X^", (b) X4//, (c) 
cos , and (d) ■03" ■ The broken curve in the cos 9s» figure is the LO QCD prediction 
for qq — > qq scattering. 
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Figure 16: Predicted distributions of the mass fractions described in the text for five- 
jet events produced at the Fermilab Proton-Antiproton Collider that satisfy the re- 
quirement msj > 750 GeV/c^. HERWIG predictions (points) compared with NJETS 
predictions (histograms) . 
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Figure 17: The predicted distributions of the variables describing the (74'5')-system 
for five-jet events produced at the Fermilab Proton-Antiproton CoUider that satisfy 
the requirement msj > 750 GeV/c^. HERWIG predictions (points) are compared with 
NJETS predictions (histograms) and the phase-space model predictions (curves) for 
(a) /a', (b) /fi/, (c) Xa', and (d) V^^/^,. 
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Figure 18: The predicted distributions of the variables describing the (CD)-system for 
five-jet events produced at the Fermilab Proton-Antiproton CoUider that satisfy the 
requirement m^j > 750 GeV/c^. HERWIG predictions (points) are compared with 
NJETS predictions (histograms) and the phase-space model predictions (curves) for 
(a) fc, (b) fo, (c) Xc, and (d) ^'Sn- 
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